home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
IRIX 6.2 Development Libraries
/
SGI IRIX 6.2 Development Libraries.iso
/
dist
/
complib.idb
/
usr
/
share
/
catman
/
p_man
/
cat3
/
complib
/
IMTQL2.z
/
IMTQL2
Wrap
Text File
|
1996-03-14
|
3KB
|
67 lines
____IIIIMMMMTTTTQQQQLLLL2222((((3333FFFF)))) ____IIIIMMMMTTTTQQQQLLLL2222((((3333FFFF))))
NNNNAAAAMMMMEEEE
IMTQL2, SIMTQL2 - EISPACK routine. This subroutine finds the
eigenvalues and eigenvectors of a SYMMETRIC TRIDIAGONAL matrix by the
implicit QL method. The eigenvectors of a FULL SYMMETRIC matrix can also
be found if TRED2 has been used to reduce this full matrix to
tridiagonal form.
SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee iiiimmmmttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssiiiimmmmttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
rrrreeeeaaaallll dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
On INPUT
NNNNMMMM must be set to the row dimension of two-dimensional array parameters
as declared in the calling program dimension statement.
NNNN is the order of the matrix.
DDDD contains the diagonal elements of the input matrix.
EEEE contains the subdiagonal elements of the input matrix in its last N-1
positions. E(1) is arbitrary.
ZZZZ contains the transformation matrix produced in the reduction by TRED2,
if performed. If the eigenvectors of the tridiagonal matrix are desired,
Z must contain the identity matrix. On OUTPUT
DDDD contains the eigenvalues in ASCENDING order. If an error exit is made,
the eigenvalues are correct but UNORDERED for indices 1,2,...,IERR-1.
EEEE has been destroyed.
ZZZZ contains orthonormal eigenvectors of the symmetric tridiagonal (or
full) matrix. If an error exit is made, Z contains the eigenvectors
associated with the stored eigenvalues.
IIIIEEEERRRRRRRR is set to ZERO for normal return, J if the J-th
eigenvalue has not been
determined after 30 iterations. Calls PYTHAG(A,B) for sqrt(A**2 +
B**2). Questions and comments should be directed to B. S. Garbow,
APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
PPPPaaaaggggeeee 1111